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ABSTRACT 


Orthotropic plate theory has been increasingly used to 
model plate-stiffener combinations typical of those which 
are used in ship hulls. The behaviour of a thin plate with 
large deflections is described by two nonlinear partial 
differential equations of equilibrium and compatibility. 
The orthotropic form of these equations is derived and 
solved for a large aspect ratio plate with long edges simply 
Supported and short edges fixed as boundary conditions. 
Buckling and post-buckling regions are investigated under 
combined loadings of lateral pressure, inplane edge com- 
pression and edge shear. Results are presented for virtual 
aspect ratios 1/1.5 to 1/6 and both isotropic and ortho- 
tropic plate properties in the form of design and behaviour 
Sharcts. 
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NOMENCLATURE 


y= ordinate in long direction 


x 


a 


M 
xy 


mn 


ordinate in short direction 

plate length in x-direction 

plate breadth in y-direction 

> = aspect ratio 

displacements of a point in x-, y-, and 2- 


directions respectively 


middle-plane strains 


flesuvalerigidvey GfLsecrthotroplc Pplatewine.<— 
or y-directions respectively 

effective torsional rigidity of Orenotropic 
plate 

POISSON'S Patliower OFENoOeropre place sin x= On 
y-directions respectively 

modulus of elasticity in x- or y-directions 
respectively 

modulus of elasticity in shear 


middle-plane loads per unit length 


bending moment in orthotropic plate acting 
around a line perpendicular to x- or y-axis 
respectively, per unit width 


twisting moment in orthotropic plate 


nondimensional deflection coefficient 
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xY 


2D 
xY 


mM Mik 


ea 
Airy's stress function 
HE 


hE 
4 


plates thickness mouerhiotrepic 


xy x 
Ps +. PN se + AC 
Gh? 
12 


total loads in x~ and y-directions respectively 


uniform lateral load, per unit area of plate 


constant inplane shear load, per unit length 











D 
= = = virtual aSpect ratio of orthotropic plate 
Xx 
Dy. 
x _= torsion coefficient of orthotropic plate 
{DD 
x yY 
J 
xY 
J J 
x Y 
Na? 
; = nondimensional in-plane load x-direction 
sity B) 
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Sx 
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exe EY LXxy 
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= nondimensional in-plane load 


y-direction 


pare 
Sve = nondimensional lateral load 

4 
a ‘'hD 

y 
Sab rin ; ; 
= ncondimensional edge shear load 

1D, 


effective width, x- or y-directions 
respectively 

membrane stresseS in x- and y- 
directions respectively; nondimensional 
shear stress in xy plane; nondimen- 
Sional 

bending stresses in x- and y- 
directions respectively; nondimensional 
twisting stress in xy plane; non- 
dimensional 


total stresses 
nondimenSional total stresses 
nondimensional total stresses 


maximum and minimum principal stresses, 
plane stress; nondimensional 
edge membrane stress in the x- or y- 


direction respectively 
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nondimensional bending moments 


per unit length 
nondimensional bending moments per 
unit length due to curvature in x- 
direction only (M¥ ) and y-direction 
O 
only (M* ) 
¥ MY 


O 


COCEELGCITENES 


Von Mises combined stress; non- 
dimensional 

moments of inertia of the stiffeners 
with effective plating in the x- and 
y-directions respectively 

moments of inertia of the effective 
plating alone in the x- and y- 
directions respectively 

spacings of the stiffeners extending 
in the x- and y-directions respec- 
tively 

equivalent thickness of the plate 
and the stiffeners (diffused) in the 
x- and y-directions respectively 


thickness of plate alone 
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Peel sence sis clemenmen se! 


Of major concern to the naval engineer is the reaction 
of ship bottom plates to a combination of inservice loadings. 
In gene:ral, the plate-stiffener combinations will be sub- 
jected to inplane compression or tension edge loads due to 
ship hogging or sagging, lateral hydrostatic loads and in- 
plane edge shear loads. With increasing compression and 
shear loads, a critical state will be reached where buckling 
will occur. Plate deflection will exceed plate thickness, 
however, the plate may still carry considerable loads. High 
slenderness ratio (b/h) plates are encountered in modern 
longitudinally framed ship designs. The yielding load may 
be considerably higher than the buckling load and investi- 
gation of the plate-stiffener residual strength in the post- 
buckling region is of great importance. 

Orthotropic plate theory has been increasingly applied 
to ship structures. With this analysis, the plate-stiffener 
combination is modeled as an equivalent flat plate with 
elastic properties that are different in two perpendicular 
directions. Schade [12, 13] and Mansour [9, 10]*, among 
others, have worked with this concept. 

The small deflection theory is useful only in the pre- 


buckling region where deflections of a plate are small in 


*Numbers in brackets designate references in Chapter 8. 





aS) 
comparison with its thickness. ‘the assumption of no deforma- 
tion in the middle plane of the plate is made in this case. 
If deflections are not small, ~ > > the strain in the middle 
plane is no longer negligible and must be considered. Von 
Karman's equations describe the behavior of an isotropic 
plate under combined lateral and inplane loads. The non- 
linear equations describe the behavior of the plate in both 
small and large deflection regions. 

Considerable application has been made of von Karman's 
equations. Levy [4, 6] investigated square and rectangular 
plates simply supported and subjected to both inplane edge 
loading and normal pressure. Coan [3] included the effects 
of small initial curvature for this case. Levy [5], using 
two fixed and two simply supported ends, solved the equations 
with edge loading only for an aspect ratio of 4, and studied 
pure shear for the simply supported square plate [7, 8]. 

The square plate with variations of the boundary conditions, 
subjected to edge loading only, was studied by Yamaki [17]. 
Payer [11] applied von Karman's equations to deep web frames 
and included uniform edge shear in addition to compressive 
edge stress and normal pressure. Shultz [15] investigated 
Wide plates having aspect ratios of 8 = 1.5 to 8 for con- 
ditions of a transversely framed ship. All of these 
analyses apply only to isotropic plates. Mansour [9, 10] 
extended von Karman's equations to consider slightly rec- 
tangular orthotropic plates under various boundary, loading 


and initial deflection conditions. 
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A large aspect ratio orthotropic plate with the short 
edges fixed and the long edges simply supported is considered 
in this paper. The plate-stiffener combination is subjected 
to loads of inplane edge compression, uniform edge shear and 
uniform pressure normal to the plate. Figure 1 illustrates 
the loadings and coordinate system used. 

The orthotropic form of von Karman's equations are 
solved similar to the theoretical analysis of Mansour [9, 10]. 
The IBM System 370 model 165 is employed to produce the 
numerical results. "Design" charts of deflection, effective 
width and bending moment are given for orthotropic plate 
virtual aspect ratios of 1/1.5 to 1/6 with inplane edge 
compresSive and lateral loadings. In addition, "behaviour" 
plots of plate centerline deflection and total principal 
Stresses are given for virtual aspect ratios of 1/1.5 and 
1/1.25 for various combinations of inplane edge compression, 


edge shear and lateral loadings. 





7 
= 
_ 
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2. Formulation of the Problem 


2.1 Orthotropic Material Properties 


An orthotropic material has different elastic proper- 
ties in two perpendicular directions. For plane stress in 


the xy plane, the stress-Strain relations are 


E 
> 
Ce ee ee 
x i} Vs Vy x yy 
By 
GO = ~——l(le. + vie) (1) 
y uF Vy¥y y x x 
Ty = Yuy? 


from energy symetry 


EV BV 
Sc ay. Vax 


which gives 4 independent elastic constants. 


2.2 The Rectangular Orthotropic Plate with Large Deflection 
(Figure 1) 


For thin plates with large deflection, r > 7 a satis- 
factory approximate theory makes the following assumptions: 
1. Points initially on the normal to the middle plane 

of the plate remain on the normal after bending. 


(Disregard of shear deformation.) 





1B 


2s The normal stresses tno mie divrectton st ranscveroeute 


the plate can be disregarded. 


(8) = 


: Or) 


(Plane stress where 


3. Hooke's Law relating stress-Strain applies. 


For the large deflection theory, deformation in the 


middle plane of the plate must be considered. 


These strain 


components therefore include the effect of deflection and 


are approximately 


5 (2) 2 


yo 


Equilibrium of forces and moments 


on 


plate produces the equilibrium equations 





oN, 
__ pe) a 
a oy 
oN 
ON, EST 
oY OX 
and 
97M 92M 92M 
X_ 95 XY 4 4 
9x2 oy ay? 
_ - law 9 “w 
=<-(g + i ile 2Nvaxay * 


N 
¥ 


(2) 


an element of a 


CS) 





(4) 
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Introducing the familiar Airy's stress function F which 


Satisfies equations (3) 


Roe lig oe Jaane 4 me JCF F 
N DAN > Gaeta ae Ee =) 








From substitution of large deflection strains (2) into 
the stress-strain relations (1) and using the definition of 
bending and twisting moments, the moment-curvature relations 


for an orthotropic material are derived 








2 2 
m= -p (2! + y Bw 
eae ye ay? 
2 2 
M, = -D (QB + y oY) : (6) 
Y ay 7 yx 2 
Moy 2c OX0Y 


where the rigidity coefficients are defined as 


E h? 
= x 
X 12(1-v,v,) 


Eon 
dD. z 12(1l-v_v_.) 
x Yy 
and 
Bech 
C = Ts 


Substitution of equations (6) and equations (5) into equation 


(4), the "equilibrium" equation is obtained as 
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paw 4 ap atv 4 p atv 
ay V9 x 2a y2 Yay 
Ci 
= 2 2 2 Z 2 Zz 
= gq + fie Chae 99 : mie oF . otw 
: 2 2 dX0dY dXdY 2 2 
oy OX OX oY 
where 
2D =Dy +Dy + 4C 


XY ay y= 


To obtain the "compatibility" equation, the strain equations 
(2) are first differentiated and combined to eliminate dis- 
placements u and v 
2 2 2 
o“€ 7 3 Yxy OE 


RY 2;_2°W,2 _ o2W . a¢w 
OXdY 2 axoy? (8) 








ox 


Differentiation of equations (1), substitution into equation 
(8) along with (5) and uSing the equilibrium equations (3), 


the "compatibility" equation is obtained as 





h b b 
yo + ag. —2 = + 5 dF 
yx Yax2yy2 Yay" 
Sy pak (9) 
aXdY 5x2 Jy? 
where 
ale i 
x Eh y Eh 
and 


co ~ GR tle = oe 
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Equations (7) and (9) are fourth-order non-linear 
partial-differential equations that describe both small and 
large deflections of an orthotropic plate. This includes the 
buckling and post-buckling behavior of the stiffened plate. 
Substitution of the isotropic material properties 


E 
x y Xx y 2(1+v) 





result in the familiar von Karman's eguations. 

Solution of these equations with 16 boundary conditions 
results in the two functions F and w. The boundary con- 
ditions are eight support conditions and eight edge load/ 


displacement conditions. 


2.3 Bending and Membrane Stress 


The membrane stresses are determined from equations (5) 








:- _ Ny Db s6-r 
o, ROT Ri 
x dy? 
N 2 
SY ax? 
XY h h dxoy 


Noting that the maximum normal stress acts on those sections 
parallel to the xz or yz planes, and using equations (1), 


(2) and (6) the bending and shear stresses are obtained as 
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E M 
a. — 5 ae — 
bx GR ee: 
x Y x 
E M 
0 en Ime A —- — 
by Gli vey) D., 
Mey 
{ — ew 
T oxy ZG C 
The maximum bending and shear stiesseS occuring at z = +3 
gives 
My. ny Mey 
O = +6— ; = +6— ; = 76 IAL 
bx h2 °by h2 Toxy ah h?2 ( ) 
where My and BL and Hh oy are given in equations (6). The 
total stresses are the sums 
— su! 
Ctx Ox) atx 
O =o +t (12) 


Oo 
ty y by 


txy ay bxy 


2.4 Load Boundary Conditions 


The following most general conditions state that the 
edges are subjected to an inplane average compressive load 
per unit length in the y-direction and x-direction of 
magnitudes N, and No respectively. Additionally, all edges 


are subjected to a constant inplane shear load per unit 


length of magnitude §. Referring to Figure l, 





BE, 








at x = tS 
=: O°F os 
Nyy G f Oxay 
Inplane load resultant Py is (13a) 
By 2,0 
ee ol oF ay = -ib 
=b/2 oy- 
at y = £2 
_= _0?F _ iz 
Nox =o 4 aXOY S 
Inplane load resultant ey is (13b) 
a/2 2 - 
Pp = f o’F dx = -Na 
yy -a/2 yx 2 y 


Pao cupport Conditions 


Edges simply supported at x = + 


NESS 


Deflection equals zero: w= 0 


External moments equal zero: (14a) 








9x 2 y ay 


Edges clamped at y =12 
Deflection equals zero: w = 0 
OW 
Slope equals zero: = = 0 (14b) 


dy 
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3. Theoretical Analysis 


3.1 Analysis Procedure 


The solution method used is identical to that of 
Mansour [9, 10] which is an extenSion of Levy [4], Coan [3] 


and Yamaki [17]. Briefly, the outline of the procedure is: 


1. Express the deflection of the plate satisfying the 
Support boundary conditions, in a double trigono- 
metric series choosing only a finite number of 
terms. w will be a function of unknown nondimen- 


Sional coefficients De 


2. Substitution of this expression into the "compati- 
bility" equation (9) results ina fourth order 
partial-differential equation with the Airy stress 
function F and quadratic functions of the coeffi- 


cient b expressed as coefficients C__. 
mn pq 


3. A solution for the stress function is assumed which 
Satisfies equation (9) and the load boundary con- 
ditions. F will be expressed as functions of 
unknown coefficients De which in turn are functions 


Of CoOcrtficrents Gus. 
pq 


A TO determine ene uniknown coecrr Teens ban! Galerkin's 


Method is applied to the “equilibrium" equation (7) 





Ins 


with substitution for w and F by their appropriate 
expressions. Using the orthogonality properties of 
the trigonometric functions, a set of simultaneous, 
non-linear algebraic equations involving cubic 


products of the coefficients Dine results. 


eee Solution 


The deflection of the plate surface can be expressed in 
the form 


4 a ve Ey) g,, fy) CS) 


Satisfying the boundary conditions (14), the deflection terms 


take the form 


fa (x) 


i 
Q 
0 
| 
* 
3 

I 
kK 
WW 
5) 


SOP) = eyes + cosa eno i= wa 2a See 


Differentiation of expression (15) and substitution into 


the “compatibility" equation (9) 





\ 4 i 
Jo + 29 eee ee 
ec oo ee Oa 
=h?{zp Eb (BN (287) sin®™x sin caTy }? 
Bf? Bae lite 2nt mi 
h {s x Dearie 1) + cos=— ay) (= ey cos ees 


{5 


jan | oa | 


Zr 2 — 
Sia (Cre) cos——x cosy} 
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ie 4b... beng sins—x sinStty 
og 





mn i} 





yo DD 4b_b; .m?j*cos=—x 
ro os J 


= [ie ee cos tty) cos x cosy} (16) 


ste eee ee eS re ere 


Nn, J} = 1, Z. et ane 


which can be expressed in the form 





4 l l 
x ae oe a ay. 
ax" ¥ 9x2 ay? ae 
He oe 2pT 





2qT 
—=-x cos==—y (17) 
Dyer Cd = 205 lipo2.as 


where one are quadratic functions of the nondimensional 
coefficients b_.. 
mn 


A particular solution for (17) is assumed in the form 


x > cos*PTx cossaty (18) 


2, 
p Pq bq 


Substitution of (18) into equation (17) defines coefficient 


ete as 
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exer = bol oa Dante 4 Soe 
16 (Jp gt + 25,..,P 2 Fag CR 5 ae ) 


where 


[oO 


The stress function F to satisfy the load boundary conditions 


(13) and equation (9) is 


— pre: x 
— = —_— — —= ee 2 = 


The method of B.G. Galerkin is applied to the "equilib- 
rium" equation (7) to determine the unknown coefficients ban’ 
Galerkin's condition requires that the following equation be 


satisfied by all functions £ (x) gly) 














a/2 b/2 \ \ \ 
fn be [D2 ¥ + 2, —-¥_ +p 2* 
ax Y9x29y? Yoy* 
2 2 2 2 2 2 
Ge ne oa gee, de ee (21) 





2 yx 2 dxXxdyY daxdY Jy 2 ay - 


- f(x) goly) dxdy = 0 


where £ x) ge) is given be equation (15). Equation (21) 
may also be obtained by applying the principle of virtual 
work. 

Substitution of w and F, expressions (15) and (20) 


respectively, into equation (21) and using the orthogonality 
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properties of the trigonometric functions, there results 


ny 2 
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where a is the Kronecker delta, me = 2002 © tee eee 
rf 
aes «(0 ae O if n = Ss and 
‘’n-s 
ae : m+xr 
A ar a eee ale ——— odd 
or 
eee BES saa 
mx eek 8 2 


Equation (21) is identical, except for the addition of the 
two shear and two N load terms, to that of Mansour [9] case 
(i) and, if isotropic material properties are substituted, 
reduces to Yamaki [17] case IIIb for zero initial curvature. 
Substitution for coefficients Pe in equation (22) 

results in a set of simultaneous, non-linear algebraic 
equations involving cubic products of the coefficients bon 
Solution of these equations gives the coefficients b and 


mn 
hence defines the functions w and F. 


eeoeebending and Membrane Stresses 


The membrane stresses are now determined from equations 


meu), (18) and (20). These are 


Ny. oa 2ptT 2qT 
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= _s _ pr can : pn 2qT 
Lg 3 isle? z d a —=—) (—=—) sin—x sin=—y (23) 
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The membrane stresses will be represented in the most 
general form, the nondimensional form used by Mansour [9]. 


Using equation (19) 
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The bending and twisting stresses from equation (ll) are 
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The normal bending stresses may be expressed as functions of 


nondimensional bending moments M* and ue where 
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The bending stresses as functions of the nondimensional bending 


moments are 
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iie nondimensional moment due to curvature in the x-direction 


only (Me ) and that due to curvature in the y-direction only 
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The bending stresses may be expressed as nondimensional bend- 
ing stresses in a form compatible to equations (24) and (25). 
Using equations (15) and (26), the nondimensional bending 


and twisting stresses are 
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The nondimensional total stresses are therefore, from (12) 


= o* + o% 
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3.4 Effective Width 


The effective width of a rectangular plate which has 
buckled is defined as that width of a uniformly stressed 
phantom plate of the same thickness stressed to the same 
maximum stress and sustaining the same total force as the 


real plate. The effective widths a. and Der in the y and x 
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directions respectively are hence defined as 


N 

ty, = = —H (34a) 
N,. 

Doe = b h (34b) 


where ove is the edge membrane stress in the y-direction 


(o. at x= t 5) and oO is the edge membrane stress in the 


y xe 
x-direction (o, at y = +t 3) . Substitution of (25) evaluated 
at x = t 5 into (34a) gives 
a N* 
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Clearly, the effective widths a. and bd, are not constant over 
the y and x-directions respectively. Neglecting the small 
amount of change due to the periodic terms in y or x, taking 
only the average values into account [15], the effective 


widths take the form 


a N* 
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3.5 Principal Stresses 


For the state of plane stress, the maximum and minimum 


principal stresses, o, and o are given by 


1 DN 
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Substituting the nondimensional total stresses, equations 


(33), the nondimensSional principal total stresses are 
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4. Numerical Solution 


foe, Analytic Solution 


Shultz [15] found that eight deflection terms were 
sufficient to describe a simply supported plate up to an 
aspect ratio of 8. In this solution eight deflection terms 
in the y-direction and one in the x-direction were assumed. 
From equation (15) with m = 1 and n = 1,2,3...8, the 


deflection of the plate takes the form 


(-1+cos- ty) 


bes TT 27 
w=h cos=x[b,, (1tcos+y) + b 5 
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67 81 107 
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LOA, 147 Mon 
+ b,,(-l+cos=y) + b, 4 (ltcos~—y) + Dig! 1+cos-—y) 


Similarly, the right hand side of equation (16) is ex- 
panded for m, i= 1; n,j = 1,2,3...8. The coefficients Cog 
are determined by collecting terms on the right hand side of 
equation (16) and matching the coefficients to the series 
on the right hand side of equation (17) for p = Q, 1; 
q= 0, 1, 2...16. The coefficients Me and is are listed 
in Table 1 where the coefficients Oye are given by 


eee (Ek (38) 
Set Gl q 
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Ga Oem 2 eer Lb 
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Substitution of coefficients oad into equation (22) 
meme equation (19) and with m, r = 1; n, Ss = 1,2,3....8 


produces the following equation after considerable simpli- 
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and 
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Equation (39) gives a set of eight simultaneous, non- 
linear equations for the eight coefficients bays Dios bi 3) 


b ane. ba 


b 17! 18 


b b 
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4.2 Computer Solution 


Equation (39) is programmed on the IBM System 370 model 
165 uSing Fortran IV level Gl. The subroutine ZEROIN from 
joe M.1.T. MATHLIB program library is used to solve the 
system of simultaneous, non-linear equations. This sub- 
program uses an iterative method of solution and convergence 
is achieved when the difference between two successive values 
of aan is less than Nee 

The foregoing analysis is for the completely general 
case with all the loadings as indicated in figure 1, however, 
the computer solutions were restricted to certain specific 


cases. 
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BaSically, two programs were written to calculate the 
deflection coefficients. The first program produces curves 
of deflection at the center of the plate, effective width 
and bending moment in the y-direction at the center of the 
fixed supports versus inplane edge compressive load NY with 
Q* aS a parameter, ne = 0 and S* fixed. Curves of this 
nature, for a range of orthotropic parameters y and Mn, 
result in a set of design charts. The other program permits 
investigation of the behaviour for any set of loading con- 
ditions with Ny = 0. Deflection along the centerline x = 0 
and total principal stresses along the centerlines x = 0 


and y = 0 are plotted. 
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5. Results 


eel Design Charts 

The design program was used to generate a set of de- 
sign charts with NY = S* = 0. Virtual plate aspect ratios 
of 1/1.5, 1/2, 1/4 and 1/6 are presented. For each aspect 
ratio a complete set of curves for orthotropic material 


parameters y = n = 1.0 and y = 4.0, n = 0.5 are given. 


Description of design charts: 
5.1.1 Charts of Deflection 

The nondimensional deflection at the center of the 
plate, from equation (15) with x = y = 0, is plotted versus 
the inplane compressive edge load Ne with the lateral load 
Q* as a parameter. Figures 2,3,8,9,14,15,20 and 21 give the 


deflection. 


Sel.e2 Charts of Effective Width 

The nondimensional effective width given by equation 
(36b) is plotted versus the inplane compressive edge load 
N¥ with the lateral load Q* as a parameter. Figures 4,5,10, 


11,16,17,22 and 23 give the effective width. 


Bees Charts of Bending Moment 


The nondimensional bending moment i’ at the middle of 
the supports y = + 2 and x = 0, equation (27), is plotted 
versus the inplane compressive edge load Ne with the lateral 


load Q* as a parameter. Figures 6,7,12,13,18,19,24 and 25 


give the bending moment. 





| 
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5.2 Behaviour Charts 

The behaviour program was used to investigate a few 
general sets of loading conditions with NY = 0. Virtual 
plate aspect ratios of 1/1.5 and 1/1.25 are presented for 


orthotropic material properties y = n = 1.0. 


Description of behaviour charts: 
5.2.1 Charts of Deflection 

The nondimensional deflection, from equation (15), along 
the centerline of the plate (x = 0) is plotted versus y/b 
with either S* or NY aS a parameter and Q* fixed. Figures 


26,27,36,37,46, and 47 give the deflection. 


5.2.2 Charts of Maximum and Minimum Total Principal Stresses 
The nondimensional total principal stresses, equation 
(37), along the centerlines on the bottom surface of the 
Plate (z= -2) are plotted versus y/b (x = 0) and x/a (y = 0) 
with either S* or NY aS a parameter and Q* fixed. Figures 
28-31, 38-41, and 48-51 give the stresses for x = 0 and 


Figures 32-35, 42-45, and 52-55 give the stresses for y = 0. 





6.) Discussion GLoResultS 


fe, Orthotropic Properties 


Reference [10] contains a complete discussion of the 
rigidity coefficients Dy» De and Bhs and compliance coef- 


fecLtents Jy J and Tey which are used to express the non- 


SY 
dimensional parameters and orthotropic material coeffi- 
cients. Possible approximate formulas for their calculatic 


are given in that reference. The two extreme cases of 


orthotropic plate coefficients are the isotropic case wher¢ 


mM = y = 1.0 and the grid (intersecting beams, no plate) 


n QO and y = ~. Most plate-stiffener combinations fall 
between n = y = 1.0 and n = 0.5, y = 4.0; hence the choice 
of parameters for the design charts. | 
Comparison of Figures 2 to 25 as to the nondimensional 
load parameters NY and Q* indicates as p decreases the 
assigned values of NY decrease and Q* increase. Nhe and Q* 
are indirect functions of the aspect ratio since they are 
nondimensionalized by the plate edges a and b respectively. 
The ranges of values used in the design curves were predi- 


cated on achieving a certain minimum value of effective 


Peeadth in the calculations. 


6.2 Comparison with Existing Solutions 


The first set of design curves to be produced were 
those for the case 9p = y = n= 1.0, Figures 61,62 and 63, 
FOr NY with Q* as a parameter and S* = N* = 0. This choice 


of parameters permitted comparison with the solution obtaince 
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by Mansour [10]. The number of coefficients retained in 
this analysis (only one in the x-direction vice two by 
Mansour) does not favor the solution for the square plate 
case, however, the results of Figures 61, 62 and 63 agree 
Beewithin 3% of the figures 1, 2 and 3 of Mansour [10]. 
The deflection coefficients obtained in reference [10] show 
that the one term assumption only neglects terms that are 
less than 5 1/2% of the primary deflection term Diy which 
accounts for the reasonable agreement in this case. 

For the loading conditions of the design curves, 
Figures 2-23, comparison was made with the solution of 
Levy [5]. Levy calculated deflection at the center of the 
plate and effective width for the isotropic plate of 
B = 1.5, Q* = S* = 0 and the same boundary conditions as 
this investigation. Comparison of Levy's results, Figure 3. 
and Table 11 of reference [5], with the curves for Q* = 0 
from Figures 2 and 3 is graphically displayed in Figures 64 


and 65. Correlation of the two results is very good. 


6.3 Design Charts 
feo... Charts of Deflection 

The deflection at the center of the plate is always the 
maximum deflection since only the first deflection mode is 
encountered in these cases. With reference to the curve 
Q* = 0 in the deflection charts, the plate remains unde- 
flected until the critical compressive load N* is reached. 


Increasing be further buckles the plate and plate deflection 





45 

increases rapidly. Hence, intersection of the curve 0* = 0 
and the NY axis gives the lowest buckling load. The non- 
linearity of the curves is evident. In the small deflection 
range, = < > a doubling of the lateral load Q* will double 
the deflection, however, in the large deflection range the 
lateral load effect becomes nonlinear with decreasing incre- 
ments of plate deflection for increasing additions of load. 

As op goes from 1/1.5 to 1/6 the magnitude of the maxi- 
mum deflection of the plate is considerably reduced. The 
plate becomes more of a beam supported all around and hence 
"stiffer" than the slightly rectangular plate. The effect 
on the deflection of increasing Q* as po > 1/6 becomes less. 

Comparison of buckling loads and the range of Ne values 
over the four aspect ratios shows that for p = 1/1.5 the 
maximum value of Nee considered is four times the buckling 
load while it is less than twice the buckling load for 
o = 1/6. The reason for this selection of Ne values is 
‘that the range where the solution is reliable decreases as 
o + 1/6. Convergence to the correct solution becomes in- 
creasingly more difficult. In addition, the 8 term solution 
begins to become insufficient for p = 1/6 where the eight 
term, Bigs has increased to a value almost 6% of the primary 


term, bias for the highest set of Ne and Q* considered. 


Oeo.2 Charts of Effective Width 
By edet titi ton,» tne plate is TUE ser Lect iy cummin ca iiayelng 


the external compressive edge load until it buckles. With 
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reference to the effective width figures, b/b = 1] for 0* & 
0 until the buckling load is reached. Extensions of the 
curves for Q* # 0 at values of NY less than the buckling 
load are meaningless and should be ignored. After buckling, 
the effectiveness drops off as the plate deflects. Compari- 
son of different aspect ratios shows that the effective 
width decreases more rapidly and takes on significantly 


lower values ap > 1/6. 


6.3.3 Charts of Bending Moment 
The bending moment - at the middle of the fixed sup- 


ports y = + D 1s calculated because it is expected to be 


2 
large. ee 1s equal to My at these points since the curva- 
O 
ture along these supports in the x-direction is zero (Me =0), 


fo) 
hence equation (28) is plotted. With reference to the 


bending moment figures, the bending moment is zero until 
the plate deflects. Increasing Q* increases the bending 
moment in all cases with the nonlinearity associated with 


large deflections. 


6.4 Behaviour Charts 

Unfortunately, there are no existing solutions to which 
the results of shear loading may be compared for this par- 
ticular set of boundary conditions. Levy [7,8] and Payer 
[11] considered the simply supported isotropic case. How- 
ever, the results obtained are consistent with those antici- 


pated in the magnitude and shape of the center line 
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deflection curve and values of principal stresses. 

Results are given Fa Se foro. =" 1/71... 5 ana We 2S 
Investigation of the solution revealed that for virtual 
aspect ratios greater than about 1/1.75, the plate will have 
one symmetric buckle about the diagonal while for smaller 
values of 0 antisymmetric buckling shapes are experienced. 
The choice of one term in the x-direction limits the solution 
to some degree in the former case*, however, it cannot begin 
to approximate the true solution in the later case. Hence, 
the narrow range of investigation for op. 

Additionally, with reference to the nondimensionalizing 
of the total stresses, sections 3.3 and 3.5, to present the 
results in a concise and meaningful form, the solution is 
restricted to the isotropic case where y = n = 1.0. This is 
a result of the isotropic aspect ratio 8 appearing expli- 
citly in the nondimensional equations. Nondimensionalizing 
the bending and membrane stresses to eliminate g will result 
ina set of incompatable nondimensional stresses which must 
be plotted and evaluated separately, thus increasing the 
number of charts required for each specific case. 

Initially, the total principal stresses at both the top 
and bottom (z = + 2 respectively) of the plate were investi- 
gated. Figures 56-60 show the results for p = 0.667, y=n=#= 
1.0, S* = 14 (above critical load) and N* = Q* = 0. Using 


the Von Mises yield criteria for the biaxial case, 0, = OF 


*Note: The assumption of the one team solution in the x- 
direction will model a symmetric buckle about the plate 
centerline, not the diagonal. 
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(ox) * = (oF) * + ot) -~ ot go (40) 


* * 
WZ 
a combined stress may be calculated to evaluate the relative 
importance of the total stress field at the top or bottom 
of the plate. With reference to Figures 57-60, the table 


one (ox) is as follows: 





oa ee ib 
eet y= is, x=0 x= is, y=0 
Ales 12.4 Liesl 10.8 
Bottom 12.2 IS) Bee Ce 3 


The combined stresses are slightly greater at the top 
than at the bottom for the center of the plate, however, the 
magnitude and difference is considerably greater at the 
fixed supports with the bottom stresses being of higher 
value. By the nature of the upwards (+z) buckling mode, 
this is the expected result. Due to the higher combined 
stress at y = + >, the behaviour charts were plotted only 
for the principal stresses at the bottom of the plate 
(Z = - 2) . 

Figures 26-45 show the effect of increasing shear loads 


above the critical load for 9 = 1/1.5 and 1/1.25. Inplane 


edge loads NY = N* O and cases are presented for both Q* = 


yy, 

O and a finite value. Increasing S* and/or Q* increases the 
deflection and total principal stresses throughout the plate. 
Figures 46-55 show the effects of completely general 

sets of loading conditions with increasing inplane edge 


compressive loads Ng above the critical load. Comparison 


of the charts for the sub-critical (S* = 6) and super-critical 


AQ 


(S* = 14) shear loads clearly indicate the influence of 
shear on the shape of the deflection curve and values of 
principal stress. Increased shear increases deflection and 
principal stresses. In addition, the highest maximum 
principal stresses are observed to occur at the fixed 


Support, y= +t >. 


6.5 Examples Demonstrating Use of the Charts 


The following examples are given to demonstrate the use 
of the design and behaviour charts. The approximate formu- 
las for calculation of the rigidity and compliance coeffi- 
cients are discussed in detail in reference [10]. The 


formulas are listed in Table 2. 


6.5.1 Design Example 


Consider the following orthotropic characteristics and 
nondimensional loads for the stiffened plate: 


o- O.o, y = 4.0, n= 0.5, NY = 4.0, Q* = S* = 0 


From Figure 9 





The critical load N. = 1.8 
as a* 


The center of the plate deflection 
w = 1.64h 
where h is the average stiffener's depth plus the 
plate thickness. 
From Figure 11 


The effective width b, = 0.505b 
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Edge membrane stress (equat:ions 24,34b and 35b) is 


N* 
De 


b - it 4.0 
oi a) Weare sao | ORS OS EOIAS 
e e 12(1-v2) : : 

2 
0 730 
e a’ J 
y 


Peom Figure 13 
The bending moment at the middle of the fixed support 
M* = -102 
y 


The maximum bending stress at this point (equation 29) 





is 
6M* hD D 
o, = + = 62 
y aye Jee let 


6.5.2 Behaviour Example 


Consider the following orthotropic characteristics and 
nondimensional loads for an isotropic plate: 


o0 = 0.667, y=ne=1.0, N* =12.077 0°: — 6.0, S* —==6.0 


From Figure 46 
The maximum deflection at the center of the plate w = 
74 hh 

From Figures 48 and 50 


The maximum principal stresses occur at the fixed sup- 
b 


ports y = + zr X= 0 and are 
2 
1, = 6.590 fb 
Deu 
x 
2 
oe = 65a 
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and the combined stress 


2 
= 54.6 AUnSLSTER 
bow 
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7. Conclusions and Recommendations 


The problem of a large aspect ratio orthotropic plate 
with the short edges fixed and long edges simply supported 
has been modeled for the case of inplane edge compressive 
and lateral loads and subcritical edge shear loads. The 
design charts provide valuable information of use in estéb- 
lishing design criteria for the plate-stiffener combination. 

From the behaviour Figures 46 to 55, it is evident that 
the highest principal stresses occur at the fixed supports 
y= dt >. The design Figures 2 to 25 do not present informa- 
tion to permit the calculation of the membrane stresses at 
the fixed supports. Design charts of oO, , edge membrane 
stress in the y-direction (denominator of equation (36a)), 
should be produced. 

As noted in the results, the one term assumption in the 
x-direction will not model the deflection of a plate with 
edge shear as the predominate loading. As the major point 
of this investigation was to model a large aspect ratio 
plate, the solution chosen was an inevitable consequence 
when the large range of p is considered. It is recommended 
that to model a large aspect ratio plate subjected to shear 
loading, close attention must be paid to the important 
deflection terms that should be included in the solution. 
Levy [8] confined himself to a 8 = 1.5 and used 14 selected 


terms in the deflection equation. 
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Additionally, the eight term solution is just suffi- 
cient to model the case for the orthotropic plate aspect 


faatero of 1/6. 
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APPENDIX A 


Details of Derivation of Basic Equations of 


Large Deflections 





1.03 
Definition of External Forces and Moments 


Shearing force per unit length parallel to x axis 


h/2 
OF =o) = az 
-h/2 *4 


Shearing force per unit length parallel to y axis 


h/2 


| a 302 
ee, 


Bending moments per unit length of sections of a plate per- 


pendicular to x and y axis respectively 


h/2 
Le = f ZO, dz 
—h/2 
hy 2 
M = ff zo az 
y ~h/2 Y 


Twisting moment per unit length of section of a plate per- 


pendicular to x and y axis respectively 


hyz 
=-f zt az 
*Y -h/2 *Y 
h/2 
M =) ee Cee 
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Strain - Displacement Relations 


Consideration of strain in the middle of the plate 





For large deflections, the effect of deflection on the 
strains has to be included in the strain-displacement re- 
lations. The accompanying figure shows a typical elongation 
of a linear element of length dx due to displacement in the 


xX and z directions from bending. 


vom Seen cle 1 (aw) 2 
Elongation = dl = 0% + ne obs 


The strain in the x direction is therefore 


~ OU 4, L(awy2 
ey = ox 7 3 5x) (a) 


The strain in the y direction is likewise 


. av 4 Law): 
ee oie (b) 


he! 


The shearing strain due to displacements u and v is 


ou OV 
r) ss ox 


he 


It can be shown that the shearing strain due to displacement 


w is 





0) S 


The total shear strain in the middle plane is 


ou OV OW qW 
oy Oe EE EE aE Cc 
Yxy * gy ~ 9x ° 9x ” DY wo 


For small angles the displacements u and v can be represerted 


by the change in displacement w as 


_ ow 

Us 2 Ox 
aw 

Nee ee 
oY 


Substitution into (a), (b) and (c) gives 








2 
Ape 
9 2w aw 
2 + =(£—) 2 
Ey at a ay? 
9 2W gw QW 
See 2? Soe, So 
Yxy ae, gx oY 


*Numbered equations correspond to those used in the main 
een t . 
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equal ial shah win) (26 Forces and Moments 


k————__——- a —______--+| 
i Se 
b 
ry: 
jee Lateral load g only 
Mi. 
M 
: Te 
| 
M OM 
x 
= ¢ M + a xX 
Z xY aX 





ane 


x i eee a 
By * Gyty 


The above diagrams are superimposed to make one loading 


condition. For equilibrium in the z-direction: 
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LF = 0 
Z 

oe dQ _ 
——daxdy + —Yaydx + gq dxdy = 0 
OX oy 

dQ 00 

ee ee 

OX i oy q (a) 


Moment equilibrium about the x-axis 


=M = 0 
SO 
aM, aM 00 7 
——*axdy _ ay dydx + Qdxdy + ay dydxdy + qdxdyO(dy) = 0 


Neglecting the last two terms since they are of higher order 


gives 
JM, oM 
a ay + Q. = 0 (b) 


Moment equilibrium about the y-axis and similarly neglecting 


small terms 





+ 
a 
Re 
! 
10 
i} 


Beeorces in the middie plane of the plate (Membrane Stresses) 
The forces acting in the middle of the plate may have a 
considerable effect on the bending of the plate and must be 


considered. 
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By 


Assuming no body forces or tangential forces 


LF = 0 
x 
oN, oN ™ 
5x dxdy + —idydx = Q 
LF = Q 
Y 


ON oN 
ay dyax a + dxdy = 0 


and from symmetry noting See = Ns 





Oe 


oN ON . 

eR yee ip 

ox oy 

3N 9N (3) 
ce See eee Aen 

oy ox 


Equations (3) are independent of equations (a), (b) and (c) 
and can be treated separately. 

Projecting normal forces N. onto the z axis and taking 
into account the bending of the plate and resulting small 
angles 


oN , 
-N ayo™ 4 (N. + se Oy (a: 2 Wax) dy 
os 1 EPs x ox ox yx 2 


Similarly projecting normal forces Se onto the z axis 


oN 2 
-N dxo@ 4 (Cees Sey one + Way) dx 
yY oy ey oY oY yy 2 


Neglecting higher than second order terms gives 


9 2W oN 
N. ——dxdy + -— ° ~dxdy (d) 
Xo 2 OX Ox 
x 
9 *w oN OW 
N.——dydx + —L - ~—dydx (e) 
Yay oy Oy 


Treating the projections of the shearing forces onto the 


Z axls the same way 


for. for Ney with slope of deflection surface in the y 


direction on the two opposite sides of the element as 


Ow OW 9 Ww 


de0) 


ON 
3 *w xy , Ow 
‘gyesey 76) eal 


op for Ne with slope of deflection surface in the x 


dixection on the two opposite sides of the element as 








s¥ and ou 4 oY ay 
Ny oaetlvax + oO sWayax 
combining (1) and (2) gives 
2n,, 3 axay + ot sNdady + — SWaxay (f) 


Differentiation of (b) with respect to y and of (c) with 
respect to x and substitution into (a) eliminates shearing 


forces on and oy 








97M 2M 2M 2M 
et ere nS eB 
9x2 OoXdY yy? OoxXdY 
Noting that Bho. = = Boy! adding expressions (da), (e) and (f) 


to the load gdxdy originally defined and making use of (3) 











2 2 2 
) M,. “ar ) Mey A ) M 
poet SES Ove 
Te yy a ey, Bs (4) 
xX XYdXdY Y 5 2 


Ape y 


BE ESL 


Moment - Curvature Relations 


Substitution of equations (2) for EL ey and ie in equations 


(1) gives* 
ox ~ I= 
oe ae 


II 


T 
xy 


eng 
-2zG——— + G 








E 2 2 

fan SM FSM)? = ay SMe Sy (S27 

ey ox? Y3y? Nas 

E ee 9 ew a ub OW) 2 5 ew fe i, (5 (a) 
Vow sy? ao Ko 2 23 ox 


d ow | OW 
OXO0Y ox OY 


Taking the moments as defined before and integrating 


h/2 

My, = f ZO az 
eee 
hye 

M = f ZO daz 


x =h/2 Y 


h/2 
=-f zt dz 


M 
yi Sayop PL 


Even functions will drop out over the integration interval 


leaving the 
for a plate 
tions. The 


*Numbers in 
TeX. 


Same expressions for moments Mie M, and My as 
undergoing pure bending and only small deflec- 
result is 


parenthesis correspond to equations in the main 


ide 








_ 9 ew ow 
M. = -D ({( + y ——) 
ci jace Yay? 
2 2 
M. = -p, (28 + y 2%) (6) 
ay7 ax? 
jo w 
M = 2C 
ay oxo Y 


where the rigicity coefficients D,. and Ph are defined as 


E h? 
eee ae 
>is PN) 
E_h? 
_ y 
Y L2ti-v vy! 


SB, 
| 


and 


Gn 
oS Sy 


"Equilibrium" Equation 


Substitution of equations (6) into (4) gives 











\ l \ 
p 2 4 2D. CEOS ee 
x ' Yax29y? Yoy* 
= 2 2 2 
=gqg+N gaat + 2N — + N ow 
9x? n4 Y Y9y? 
where 
2D =Dy +Dy. + 4C 





Jk 


Substitution of the stress function defined by equations (5) 














bi by | 4 
D- wa 2D. RAS oy 
ox 4 Y9x29y? Yay" 
at °F o¢w 40°F | 9?w 0°F O¢w (7) 
ay? yx 2 OXoY OX0Y yx 2 dy? 


‘Compatibility" Equation 


The forces Nyr Ny! and Ne in the middle plane of the 
plate depend on the strain due to bending aS well as the 
external forces applied in the xy plane. 

Again assuming no body forces and requiring load gq is 
perpendicular, equations (3) apply for equilibrium in the 


middle xy plane 


aN oN. 
“a “eG 
(3) 
oN oN 
Lo Eee SG 
oy ox 


The corresponding strain components are those of equations (2) 


. du, 1/dW) 2 
e. 7 ox 7 a 3x? 
= 9V , 1/9W)2 
f= oy 2 bay) (2) 
du OV OW ow 


a oy 
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Differentiating equations (2) and combining to eliminate 
mand v results in 


2 2 2 
Oe aly ne (a?w_y2 _ 3?w . 3?w (8) 
xy? OX dy 2 OXdY 2 dy 2 








Ox ox 


Using Hooke's Law to relate strain and forces N 


N N 
eae 
Xx hE y hE 
x yy 
Ny Ne 
e = -_t-v —X (a) 
y oN x hE, 
Nxy 
Yxy ~ he 


Differentiation of equations (a) and substitution into 
equation (8) with use of equation (3) 
3°N 3°N 37N 


4 Bed. fe es 
Xo 2 XY odXdY ye 


x 








dy 











ve 0 *w 2 9 2w 9 7*w 
(sxay 2 
AVSoC dy 
where 
1 1 
D48 Be Y Eh 
and 
Zu = Viel 


Sy, Gh x Y Vy x 
Substitution of the stress function defined by equations (5) 
4 4 
u Pa 2 ee Bs 


4 2 2 2 
9x4 Yox2ay2 9 Yay? 


J 
RO 9x7 dy? 








i) 


APPENDIX B 
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Iss 


Effective Width 





Figure b 


Consider the simple case of a rectangular plate simply 
Supported on all edges where the loaded edge remains straight 
in the plane of the plate at all times. If the load is below 
the buckling load, the stresses will be distributed evenly 
as shown in figure a. For post-buckling loads the center 
of the plate will exhibit less compressive strain than the 
edges because of large deflections of the center of the 
plate. The stress distribution is shown in figure b. The 
effective width relates the maximum stress ones? uniformly 
distributed along a phantom plate sustaining the same total 
load as the real plate, thicknesses being the same. 

Therefore 


a oO = Tad 
e max 


Where o is the average edge stress. For this particular case 





ele 


=i ah 


Qt 


and = O 


Omax ye 


Where ne is the edge membrane stress in the y-direction 


— at x = + =). Substitution into the above relation gives 


the result 


N 
ao = @a = (34a) 


N 
b.o,, = b — (34b) 
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APPENDIX 2c 


Details Of “colution 











Ar 


me Determination of CoefEicients cog 


Substitute deflection expression (37) into the right 
hand side of equation (16). The right hand side of equation 


(16) becomes* 











Term #1 
meh 7 . oT . 20 . 40 . 67 
= Seas ax lb,sinyy + 2b,Sinvy + 3b, sin+- 
es Th . LOT . 127 
+ 4b, sins y + 9b, sin-—-y + 6b -sin=y 
+ 7b sinttTy + 8b sintoty) } 
Z b 8 b 
Term #2 
= CE cos? 2x ih +. iD costly - b, + b cosy + b 
ey 2 a 1 I b 2 2 b 3 
+ b,coss@y = by + b ,cossty + be + b.cos=y - De 
+ b,cos=-"y + ba + b.cos="y = De + b,cos=-"y] + [b, coss2y 
47 67 81 107 
+4b,cosyy + 9b,cos+—-y + l6b, cosy y + 25b,.cos—— 
+36b,cos=="y + 49b. cosy ~ 64b,cos=—"y] } (a) 
*Note: since m = 1 for all cases, for convenience by . is 
represented by bi. F 
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Each of the two terms of equation (a) is multiplied out and 


like terms are collected in the form cossaty Lor 


ge— 0,1,2...16 using the function product relations 


4sinasing 2NCOS (apiece COs (ot on 


4cosacosg = 2 cos(a~-8) + 2cos(at®8) 


The first and second terms produce coefficients aXe and 2Y 


respectively. Using the function relationships 





sin’=x = > cos ix 
SO = = + scossix 
there results 
First term: 
X,cos<gty = x cosSply cos=2x 
Second term: 
y,cos “ity + y,cos=aty cos"x 


Comparison with the right hand side of equation (17), the 


coefficients C are 
Pg 


Q 
il 
Ke 
+ 
ms 
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or 


C= ¥. + Gi ex 
Sst Ge es 


Peeeepplication of the Method of B.' G. Galerkin to Equation 
(7) 


Method of B.G. Galerkin (from reference [18]) 


It is required to determine the solution of the equation 
L(w) = 0 
where L is a differential operator in two variables, which 


solution satisfies homogeneous boundary conditions. 


The approximate solution sought is in the form 


Zz i 
w(x,y) = 2 bo F (sry) 
n=1 


where FL (x+y) (n = 1,2...1) is a system of functions chosen 
to satisfy the boundary conditions and Db are undetermined 
coefficients. Consider the functions F (sry) to be linearly 
independent over a given region (no one of the functions can 
be expressed as a linear combination of the others). For 
w(x,y) to be the exact solution of the given equation, it is 
necessary that 


L(w) = 0 


1 


If L(w) is continuous, this requires that the expression 
L(w) is orthogonal to all functions of the system F(x ry) 
b 


fre 1,2...1). With i constants, b oreeb., i COnditirome 


ies 
of orthogonality can be satisfied. 
Requiring that the two functions are orthogonal over a 
region produces the following system of equations 
ffL(w) FE (x,y) dxdy 
R 


= Ae : bP, (x+y) ) F, Gry) dxdy = 0 


ial, Sot ee eee) 
from which the coefficients ba can be determined. These 
coefficients will define the solution w(x,y). 
The method of B.G. Galerkin can also be obtained from 
the principle of virtual work. 
Application of the method of Galerkin to equation (7) 
requires that the following equation be satisfied by all 


functions £_(x)g. (y) 














ae b/2 4 4 \ 
im oS [D, ee 2D. Es 5) Se 
0 oO ax" Y 9x7 dy? Y gy" 
hie Joe a ey Oe SE es 9 Ww) 
ay? 4x2 oxoy oxX0Yy yx 2 dy? 


X £ (x) 9. (y) dxdy = 0 (21) 
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Substitution of w and F, expressions (15) and (20) 








respectively 
a/2 b/2 

mm, 4 ep eae 2nmt mT 

5 gE BMP ROS) Pol OD) 5 a 





2nT, 2 ,MT, 2 2nmt mT, 
+ 2 Poa 2% (=p ) =e) b noes y cos—x 


anu 2n i mm = 
: 2% % (——) *b mnc°os py COS yx - 


ppg (BEy2(AaTy 2g ob p(-2)P**+ cos@Bty) 


pq pq mn 


° cos“aty cos x cos=PT, 


-NhHEE (SPT) 2 b nCOS PHY coss—x 


-h? y (SDT) 2 (2PT) 2 qh OS “Be cosahy cose x cos-PTx 


N MT, 2 Breese 2nt mT 
(a) i) + cos=—y]cos—=-x 


=r 


(227) (yp _sinaaty sin" 


f 
NO 
op 
2 
i 

8 ted 


ry (APT) (BY) (AEX) (78%) 5 Bey ete 


S 
+ 2h 5 & pq mn b b 


0 | 


a 
m 


sinsx sinSPTx}-[(-1) 5**+ cosS4y} cos ax dxdy = 0 
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The orthogonality properties of the trigonometric 
functions must be used to evaluate the integral of equation 


(21). By employing the function relations 


COSA + cOoSB = 2cos5(A + B) coss(A - B) 


sinA + sinB 2sin5(A + B) coss(A - B) 


The trigonometric products may be expressed in typical forms 


Be aan a 
if cos——x cos— x ax 
0 ll al 


1 a/2 
= q [2cos— (mtr) x + cos_(m-x) x + cos_(r-m) x] dx 


and 
a/2 mT rt 
f sin—-x cosS——x dx 
0 a a 


1 2/72 T aes a! 
= a [2sine (mtr) x + sins (m-r) x - sine (r-m) x] dx 


where integration 1S carried out for all m and r. By 
combining terms that integrate to Similar values, these 
integrals are equivalent to 


1 a/2 
a [2cos=|mtr|x + 2cos+|m-r|x + 2cos4|r-m|x] dx 
a a a 


and 


1 a/2 
a [2sin[|mtr|x + 2sin=|m-r|x ~ 2sin_| r-m|x]dx 


where the second two terms in each expression integrate to 


I25 


zero if r > m orm > r respectively. Manipulation of the 
trigonometric products on the right hand side of equation 
(21) into these typical relations permits the reduction of 
the equation into the convenient form of equation (22). 
Table 3 lists all the specific integrals from equation (21) 


with their respective values. 
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APPENDIX D 


COMPUTER PROGRAMS 





2H 


The computer program to solve the simultaneous equations 
(39) is written in Fortran IV Level Gl for the IBM system 
370. The program consists of the main and three subprograms-- 
MATHLIB subroutine ZEROIN, and subroutines EVAL and CPQ. 
Table 4 lists the important computer symbols used. 

ZEROIN is a Fortran IV subprogram, from the M.I.T. 
mathematical library, which computes the solution vector 
b of a set of N simultaneous, nonlinear equations 
N 
EE y (2) = 0 using double-precision arithmetic. b is obtained 
by an iterative method beginning with estimated values of 


the solution vector and iteration is performed as 


with the vector oe = ak. (b*) where F (b*) represents the 


vector of function values at the point bX and oo represents 
an approximation to the Jacobian at be Convergence is 
tested with the expression 
| K+ 


pktl _ pKi2 < 0.5 x 10712/b 


k+l) 2 


Subroutine EVAL generates the eight simultaneous 
equations and calculates the vector of function values for 
a given solution vector b. EVAL is called repeatedly by 


subroutine ZEROIN. 
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Subroutine CPQ iS in turn called by subroutine EVAL, 
with a solution vector b, to calculate the coefficients ove 
that are used in the Simultaneous equations. 

For the starting values of the solution vector b, the 
main program uSes a well guessed set of values set in the 
program. If the lateral load is not zero (Q* # 0) (see 
note 1) or the present loading conditions did not change 
Significantly from the last set, the previously computed 
solution vector b is used for the starting values. 

If ZEROIN cannot achieve convergence in 50 iterations, 
the main program alters the initial set of starting values 
and returns them to ZEROIN. Normally convergence is obtained 
Within the first 50 iterations, however, never more that a 


second set of values is required. 


Part 1 - DeSign Charts 


The main program varies the nondimensional loads Q* 
and N* or S* for a particular set of plate parameters. 
Tables of nondimensional deflection coefficients, nondimen- 
Sional deflection at the center of the plate, nondimensional 
effective width and nondimensional bending moment (y- 
direction) at the middle of the supports (y = + b/2) are 
outputs. Additionally, CALCOMP plots of the last three 


tables are produced. 
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Part, 2) =) Behar out seogueaiis 


The main program reads the input parameters for 4 sets 
of loading conditions on a particular plate and outputs 
CALCOMP plots of nondimensional plate deflection along the 
centerline (x = 0) and nondimensional total principal 
stresses along the centerlines (x = 0 and y = 0). These 
programs are completely general and can be used with any 
set of loading conditions. 

The design programs must be used with care. With 
reference to equations (36a) and (36b), the effective width 
is only defined for oe # 0 or Ne # 0 respectively. Hence, 
variation of S* with N* = 0 will give a value of zero for 
the effective width. The design programs are written so 
that no computation or plot of effective width will be 
made for any variation of S* with N* fixed. Additionally, 
deflection at the center of the plate and the y-direction 
bending moment may not be significant parameters in this 
case. It is recommended that the "behaviour" program be 


used if variation of the shear load, S*, is desired. 
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Program Descriptions 


Main Program 1 - Design (Ny = 0) 


imeut One card per set of complete design curves set up as 


follows: 


Column 


F=10 
iP 20 
24-30 


31-40 


41-50 


pal 60 


on = 70 


Parameter 


RHO 
GAMMA 
| B28 


*NXSTAR (Enter Negative Value 
to vary) 


*SSTAR (Enter negative value 
to vary) 


RANG 


RANGQ 


*Note: Either NY or S* may be varied, but not both in the 


same run. 


Output 
1. Tables 
a. Sen 
pb. Kat y=x=0 
Ch be (only if Ne is varied) 
d. 2 


My SHEA Se ey SS 0 





Ll 


2. CALCOMP Plots 
Tables b, c and d above. 
3. Values 


Ne or S* = Q to RANG in 10 increments of RANG/20 and 


Q* = 0, RANGQ/6, RANGQ/3, RANGQ 


Termination 


A blank card as the last data card is reguired. 


Main Program 2 ~- Design (NY 


E32 


0) 


[ieput one card per set of complete design cunves Jsepsumeaa. 


follows: 


Column 


PG 
EE—20 
ZA 30 


31-40 


41-50 


S60 


oi- 70 


Parameter 


RHO 
GAMMA 
ETA 


*NYSTAR (Enter negative value 
to vary) 


*SSTAR (Enter negative value 
to vary) 


RANG 


RANGQ 


*Note: Either NY or S* may be varied, but not both in the 


same run. 


Output 


1. Tables 


eee CALCOME Pilots 


a. 


lon 


Ce 


(Only Ee NY, is varied) 


x = 0 


mee: 


3. Values 


Ny or S* = 0 to RANG in 10 increments of RANG/20 and 


Q* = 0, RANGOQ/6, RANGQ/3, RANGQ 


Termination 


A blank card as the last data card is required. 
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Main BEOGhaM  oe—  2eNav1our 


Input four cards per program run, each set up as follows: 


Column Parameter 
1-10 RHO ) 
E20 GAMMA 4must be same all 4 cards 
Zi—30 ETA 
31-40 NSTAR (see note 2) 
41-50 SSTAR 
6 0 OSTAR 
Sl—70 sabre 


*Note: If TB other than zero, stresses at top of plate will 


be computed. TB = zero defaults to bottom of plate. 


Sueput 
1. Tables 
a. Ts 
love rat ee 70 
ele oF and O% at y = x = 0 
2. CALCOMP Plots 
ae Fr at x = 0 for - = — +5 


oes 0) 


dole 
1A 

wh< 
1A 
+. 

Nie 


* * = es 
i and 05 at x QO for 


| a De 
C. OF and 03 at y ORGEor > = oct 


Mermination 


A blank card as the fifth card is required. 


hole 


18S) 
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Subroutine IW 


Two subroutine programs EVAL - one for Me = 0, the 


Senex for N* = 0. The appropriate EVAL program must be 


used with the corresponding main programs. 


Subroutine ZEROIN & CPO 


Completely general subprograms. Use with all combina- 


tions of main programs and EVAL. 


NOTE 1 


If the lateral load is equal to zero (Q* = 0) a possible 
but trivial solution is that all the coefficients ban = 0. 
Care must be exercised to avoid uSing De. = 0 for starting 
values in this case. Additionally, under certain circum- 
stances the solution readily converges to zero rather than 


the appropriate non-trivial solution. The programmer should 


be aware of this possibility. 


NOTE 2 


The behaviour program may be used with either a Ors NY 


as a loading, but not both. The main program must be used 


with the appropriate subroutine EVAL. 
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Buckling and post- 
buckling behavior of a 
large aspect ratio ortho- 
tropic plate under com- 
bined loadings. 


thesA4298 
uckling and post-buckling behavior of 





